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The statistical measures of complexity C{N) defined by Lopez-Ruiz, Mancini, and 
Calbet (LMC) and T^^/siN) according to Shiner, Davison and Landsberg (SDL) are 
calculated as functions of the number of particles N for four quantum many-body 
systems, i.e. atoms, nuclei, atomic clusters, and correlated atoms in a trap (bosons). 
The strengths of disorder a and order /? are evaluated for each system by imposing 
the condition Ta^pi^N) = C{N). The proper pair {a, 13), obtained by the above 
requirement, can serve as a Pair of Order-Disorder Indices (PODI), characterizing 
quantitatively order versus disorder in any quantum system. According to the above 
classification, we assign to bosons the complexity character of disorder, to atoms 
the character of order, while nuclei and atomic clusters are (less) disordered and lie 
between them. This criterion can be used to estimate the relative contribution of 
order and disorder to complexity for other more complicated quantum systems as 
well and even classical ones, if one is able to describe them probabilistically. We 
also address the issue, whether those systems can grow in complexity as N increases. 
Our comparative study indicates that atoms are the only quantum system out of 
four, which is ordered, with the ability to self-organize. We conjecture that this is 
an information-theoretic reason that atoms are suitable as building blocks of larger 
structures of biological interest i.e. molecules and macromolecules. This is in the 
spirit of Wheeler's it from bit quote, the project to present everything about the 
Universe in terms of information theory. 



PACS numbers: 05.65. +b 02.50. -r 



2 



I. INTRODUCTION 

The question whether physical or biological systems can organize themselves without the 
intervention of an external factor, is a hot subject in the community of scientists interested 
in complexity. A practical way to answer such a question is to use a suitable definition of 
complexity and check if this quantity increases with the number of particles A^, implying 
the ability or inability to self-organize. There are several measures of complexity in the 
iterature. One of them is the algorithmic complexity according to Kolmogorov and Chaitin 
defined as the length of the shortest possible program necessary to reproduce a given 
object. The fact that a given program is indeed the shortest one, is hard to prove. In this 
paper we employ the statistical measure of complexity C, defined by Lopez-Ruiz, Mancini, 
and Calbet (LMC) [sl, which can be calculated easilyjprovided that the information content 
of a quantum system is known from previous work In the same spirit, we use the 

measure of complexity Ta,f3 according to Shiner, Davison and Landsberg (SDL) [15| . 

We intend to calculate and compare C{N) with Ta^piN) as functions of the number 
of particles N in four quantum systems, namely atomic nuclei, atoms, atomic clusters, and 
correlated atoms in a trap-bosons. Here, we use continuous probability density distributions 
of particles of quantum systems, while in [16] we obtained {a, (3) using C{Z) and Ta,i3{Z), 
taking into account discrete probability distributions of electron configurations in atoms. 

It is noted that such complexity measures, based in a probabilistic description of a quan- 
;um system, have been calculated quantitatively (as functions of Z) for the first time in 



13|, using atoms as a test bed. A previous preliminary step was the finding that Lands- 



3erg's order parameter 17l | is an increasing function of A^ for fermions and bosonic systems 
i (see also [ISJ). 

The aim of the present work is to find specific overall values of (a,/?) for each system, 
by requiring that approximately rQ,^^(A^) = C{N). We investigate the possibility whether 
the pair [a, (3), chosen in such a way (PODI: Pair of Order-Disorder Indices) is useful for 
the classification of quantum many-body systems according to the contribution of order and 
disorder to complexity. 

Our paper is organized as follows: In Sec. [TT] we define measures of information and 
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complexity, in Sec. Illll we present our numerical results and Sec. IIVI contains our conclusions. 

II. MEASURES OF INFORMATION CONTENT AND COMPLEXITY OF A 

SYSTEM 

I — I 

Shannon's information entropy [19], is defined as 

Sr = - j P(r) lnp(r)dr (1) 
Sk = - [ n{k) \nn{k)dk (2) 



where p(r) and n(k) are normalized to one density distributions in position and momentum 
spaces respectively. 5*,-, 5*^ depend on the units used to measure r and k, while the important 
quantity is the sum S = Sr + Sk, which is scale invariant i.e. independent of units. There is 
a delicate balance between conjugate spaces leading to the following results: 

First, the entropic uncertainty relation (EUR) holds of the form S = Sr + Sk > 3(l + ln7r) 
for a 3-dimensional quantum system The lower bound is attained for a Gaussian 

distribution. The above inequality is stronger than the Heisenberg uncertainty relation, 
because the right-hand side of EUR does not depend on the state, while Heisenberg's does 
depend. Also Heisenberg's inequality can be derived from EUR, while the inverse is not 
true. S represents the information content of the system in bits, if one uses logarithm with 
base 2 and in nats (natural unit of information), if the logarithm is natural. 

Second, the universal property S = a + 6 In was proposed and verified for various 

n n 

quantum systems (nuclei, atoms, atomic clusters and correlated atoms in a trap) [IQl, [Hi. 
The parameters a, b depend on the system under consideration. That property holds for 
systems of various sizes, obeying different statistics (fermions or bosons) with different num- 
bers of particles and various interactions. Specifically, the sizes of those systems range from 
the order of fermis (10^^^ cm) in nuclei, to 10^ A (10~^ cm) for bosonic systems, while the 
number of particles starts from a few to a hundred and goes up to millions (bosons). 
The Shiner, Davison and Landsberg (SDL) measure complexity ra,j3 is defined as: 

r„,;3 = A"•^]'^ = A"(i-A)^ (3) 

where 

A = and n = l-A (4) 



are the normahzed measures of disorder and order respectively, according to Landsberg [17[ , 
obeying the inequahties 0<A<l,0<r2<l. 

The parameter a represents the strength of disorder , while (5 the strength of order. We 
have a complexity measure '^a.ji of category I, if /5 = and a > 0, where complexity is an 
increasing function of disorder. In category II (a > 0, /5 > 0) '^ol,^ is a convex function, while 
in category III (a = 0, /3 > 0) F^^^ is a decreasing function of disorder A. 

The LMC measure C [3] is defined as 

C = S-D (5) 

where 5* denotes the information content stored in the system (in our case Shannon's infor- 
mation entropy sum S = Sr ^ S}^ [is]]) and D is the disequilibrium of the system. 

For a discrete probability distribution {pi}, the disequilibrium D can be defined as the 

quadratic distance of {pj} to the equiprobable (uniform) state pj = — ,i = l,2,.... Thus 

n 



Pi 

n 



i=l 

while for a continuous probability distribution p(r), it is extended in position-space as 

D^ = J p2(r)ci(r) 

and in momentum-space 

Dk = j n\k)dik) 

The extension from the discrete to the continuous case is justified in {o^. The proper com- 
bination of Dr and D^, to be inserted into (1) is 

D = Dr-Dk (6) 

In fact Dr has the dimension of inverse volume, while D}. of volume, leading to the dimen- 
sionless quantity D. 



III. NUMERICAL RESULTS 



We intend to calculate numerically the LMC measure of complexity C{N) = S ■ D and 
the SDL one, Fq_^(A^) = A" ■ as functions of the number of particles (electrons in 
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atoms, nucleons in nuclei, valence electrons in atomic clusters, and alkali atoms in bosonic 
traps). Our calculations are facilitated by our previous research in information entropy S. 
The key quantities of our calculations are the density distributions in position-space p(r) 
and momentum-space n(k), obtained as follows. 

For nuclei we performed microscopic mean-field calculations [lo!], using a density- 



dependent Skyrme force, in particular the SKIII 



orce 



2l| . For atoms we used the RHF 



(Roothaan-Hartree-Fock) electron distributions [2^, applied to calculations of atomic com- 
plexity 1^, [l4]. For atomic clusters, we employed the jellium model (Ekcard model) jlo| . 
231], 2J], while for bosons in a trap we carried out a numerical calculation of densities in 
both position- and momentum-spaces, solving numerically the non-linear Gross-Pitaevski 
equation [9i]. In all of the above calculations special care has been devoted to an accurate 
treatment of the Fourier transform of p(r) in order to obtain n(k) for each system. 

Next, the distributions p(r) and n(k) calculated as described above, are inserted into 
relations ([T]) and ([2]) giving Sr{N) and Sk{N) for each system, needed to find the sum S{N) = 
Sr{N) + Sk{N). Thus, C{N) = S{N) ■ D{N) can be obtained easily. However, in order to 
evaluate r„,^(A^) = A(Ar)" ■ (1 - A{N))^ , we need in addition A{N) = S{N)/ S^^^{N), 
where S'max(iV) = S 

rmax 

For discrete probability distributions the value of S'max, is attained for an equiprobable 
(uniform) density distribution i.e. Smax = InA^, where is the number of states with prob- 
abilities {pi}, i = 1,2, . . . N. S is minimum, Smin = 0, when only one probability is different 
from 0, specifically pi = 1 for a fixed i, while all the others vanish. For continuous prob- 
ability distributions, calculations for Smm and Smax are based to the following inequalities, 
introduced for the first time and verified for atoms by Gadre and collaborators [gI, [t], for 



atomic clusters and nuclei in 
^ (see also [251). 



Ill ] and bosonic systems (correlated alkali atoms in a trap) in 



5'rmin ^ ^ ^rmsiK (7) 

'S'fcmin ^ 'S'fc < Skma,x (8) 

'S'min ^ "S* < 'S'max (9) 

where S = Sr + Sk and p(r), n(k) are normalized to 1, i.e. / p{r) dr = 1 and 
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n(k) dk = 1. 

The lower and upper limits for density distributions normalized to one, are: 

^r.min=^(l + ln7r)-^ln(^r) 

^.max = ^(l + ln7r) + ^ln(^(r2)) (10) 



and 



^.min = ^(l + ln7r)-^ln(^(r^)) 

^™ = ^(l + ln7r) + ^ln(^T) (11) 



Thus 



= 3 (1 + Invr) - ^ In ^(r')^) (12) 

^max = 3 (1 + Invr) + ^ In (^(r')r^ (13) 

where (r^) is the mean square radius of the system and T its kinetic energy. We use, as an 
input in A = -5^, the values of S'max according to relation ffT3l) . 

'-'max ^ ^ I 1 I 1 

We repeated and verified the above inequalities for atoms in [l3|, [ij] and next, we 
checked their validity numerically for nuclei, atomic clusters and bosonic traps in 9|, 10|. 
We note that (r^) and T, can be calculated employing p(r) and n(k). Thus, we obtain the 
function S'max (A^) for each quantum system. 

Following the procedure described above, we calculate and plot the functions S{N), 
Sraax{N), C{N) and Ta,i3{N) (for various values {a, (3)), and D{N), for all four quantum 
systems under consideration. Our results are presented for atoms in Fig. [H atomic clusters 
in Fig. [21 nuclei in Fig. [3] and bosons in a trap in Fig. HI We note that in Fig. [T] we 
present both the dependence of C on all values of N and separately the dependence of C 
only on values of N for closed shells atoms. The data of the latter figure is used in the 
fitting of Ta^/3{N) to C{N). We also plot linear fitted expressions for S and S'max of the form 
S' = a + 6 In A^. In order to choose the specific values (a, /3), we employ a new prescription, 
described below. 

We begin with the requirement that the curves Ta^jsiN) and C{N) should coincide (ap- 
proximately) or show the same pattern, for a proper pair of values (a,/3). We quantify the 
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similarity between the two curves with a norm (Q — (rQ,_^)j)^, which should be zero in 

i=l 

the ideal case, where the two curves are exactly the same. A comparison of rQ,^^(A^) with 
C{N) has been suggested for the first time in [1^], where we obtained by inspection of the 
figures that roughly a ~ and /5 ~ 4 (a < /?) and has been observed that complexity 
increases with N = Z, based on the trend of closed shells. 

We limit our search for the optimal pair {a,f3) only for the regions of the a — f3 plane, 
where the behavior of the complexity measure Fq,^^ is the same with the corresponding 
behavior of C{N). Thus, if for example, C{N) is a decreasing function of A^, then we search 
for the proper pair {a, (3) only in the region where Fq,^^ is decreasing with A^. The trend of 
complexity is dictated by C{N). 

The resulting pair of values is named Pair of Order- Disorder Indices (PODI). We 

hope that it will serve as a useful tool in order to classify quantum (or classical) many body 
systems, according to the contribution of order versus disorder contributions to complexity. 

In Fig. [5l we present our results for all four quantum systems under consideration in the 
a — (3 plane with (0 < a < 20, and < /? < 20), with mesh points 10^^ for a and (3. A 
decreasing behavior is denoted by white color, a convex behavior by grey and an increasing 
trend by black. 

Qualitatively, the structure of the corresponding plots m. a — (3 plane and the shape of 
the complexity regions are the same for all considered quantum systems. A monotonically 
decreasing region is followed by a convex one and then we have a monotonically increasing 
area. 

Our results are summarized in Table [H where in addition we present the approximate 
linear limits of the three respective complexity regions (decreasing, convex, increasing). 
Of all four quantum systems, atoms and bosons are exceptional, in the following sense. 

• Atoms: In atoms there are very hard oscillations of complexity measures with A^ (= Z, 
atomic number). This is displayed in the corresponding (fourth) plot of Fig. [1], where 
we employ, just for the sake of comparison, the values a = and (3 = 0.42, originating 
from the fitting C{N) = Ta^/jlN) for closed shells (fifth plot of Fig. [1]). It is seen 
that the overall behavior can not be characterized with certainty as monotonically 
increasing, decreasing or convex. It is clear that our proposed two-step procedure to 
find [a, (3) and the complexity character i.e. first to observe the general trend of C{N) 
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Quantum 
System 


Trend of Complexity C{N) 
Decreasing Convex Increasing 


PODI pair 

a /3 


Complexity 
Character 


Atoms 
(norm=4.095) 


a > 9.09 p a< 9.09 /3 a< 3.33 /? 
& a > 3.33/3 


0.42 


a < /3, a = 
Order 


Atomic Clusters 
(norm=2.88 x 10"^) 


a > 2.90/3 a < 2.90/3 a < 1.27/3 
& a > 1.27/3 


10.57 1.17 


a> f3 
Disorder 


Nuclei 
(norm=2.48 x IQ-^) 


a > 2.47/3 a < 2.47/3 a < 1.25/3 
& a > 1.25/3 


10.98 


a > /3,/3 = 
Disorder 


Bosons 
(norm=5.6 x lO^'^) 


a > 250 /3 a < 250 /3 a < 15.39 /3 
& a > 15.39/3 


120 0.44 


a » /3 
Disorder 



TABLE I: Analysis of results for the PODI pair {a, (3) and complexity character of various 

quantum systems. 



and then fit Ta^jsiN) to C{N) is not straightforward for atoms, if we insist to employ 
the dependence of C on all the values of the atomic number N = Z. Is seems that 
details of our method should be adjusted, to some extent, to available data and/or the 
specific system under consideration. 

On the other hand, one can find the complexity character based on specific character- 
istics of the system that may simplify the whole picture (something that can be useful 
for further studies of complexity especially in classical chaotic systems). In the case 
of the atoms, such a simplification should be the choice of the closed shells atoms. 
Thus, we organize our results to define decreasing, convex and increasing behavior 
using the trend of the closed shells atoms {Z = 2, 10, 18, 36, 54). This choice for atoms 
implies an increasing trend for C{N) and a optimal pair (a,/3) = (0,0.42). The same 



;rend has 

isl: 



14 



3een observed in a series of papers employing various models and methods 



16 



30 



3l|. 



Bosons: The increasing area dominates, while a very narrow decreasing area exists for 
very small values of (3. The optimal pair is (20,0.08). Extending our calculations in 
the a — (3 plane for < a < 120, and < /3 < 120 (we do not display this extension 
in Fig. SD, we see that the best SDL-LMC fit is attained for a = 120 and (3 = 0.44. 
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Regardless of the exact proper value of a, which is certainly larger than 120, we are 
sure that a ^ (3. This is a clear indication that a bosonic system is to a high degree 
disordered, as expected from the absence of the Pauli principle. 

We observe in Fig. [IM] and the norms of Table [H that the curves C{N) and Va,p{N) 
almost coincide for the optimal in clusters and nuclei, while for atoms and bosons, 

they show a similar pattern. This is the best result we can obtain. 

It is concluded that three out of the four considered quantum systems (atomic clusters, 
nuclei and bosons) cannot grow in complexity or organize themselves as the number of 
particles increases, because C{N) = Ta^jsiN) is a decreasing function of A^. The case of 
atoms is different, as discussed above. The question is open whether those systems can show 
self-organization (organized complexity), when they form more complicated structures. At 
the moment, our indicative result is that for atoms complexity increases with the atomic 
number. This is interesting and promising for future research, because nature employs atoms 
as the basic building blocks of molecules and macromolecules related to biology. The next 
step is obviously to change their environment by influencing them by an external factor. In 
12] it was found that increasing the magnetic field applied on excited electrons of atoms, the 
information content of electrons in specific states is exchanged, while there is simultaneously 
an energy level avoided crossing of the corresponding states. An other idea is to study the 



effect of confinement 2m\, simulating, in a way, the effect of the environment. 

The most obvious objection to SDL and LMC measures of complexity Ta^/B and C respec- 
tively, is that any function of 5* e.g. Ta^f^i^), (A = S'/5'max) does not give anything new, 
compared with S. It is still a function of S. This might be true for some cases, but in our 
specific treatment of quantum many-body systems with A^ particles, the function S'(A^) is 
different than S'max(^), providing a non-trivial dependence of Ta^i3{A{N)) on A^. 



IV. CONCLUSIONS 

It is seen from Table H] that for bosonic systems a is very large and (3 very small, (3, 
which implies that they are disordered, as expected intuitively, because of the absence of 
the Pauli exclusion principle. Atoms can be considered as ordered {a = 0). This seems 
reasonable and might be expected for a system created by filling with electrons, one by 
one, well defined orbitals in the Coulomb field of the nucleus. Complexity measures are 
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oscillating hard in atoms, but our experience leads us to the realization that the proper way 
to assess complexity is to define the overall trend by observing the closed shells atoms. A 
similar conclusion has been drawn for atoms in [14] and recently in |l^, where a = 0.085 and 
f] = 1.015, obtained with discrete (fractional) occupation probabilities of atomic orbitals, 
indicating an ordered system. In the same spirit, nuclei and atomic clusters are (less) 
disordered systems {a > (3) lying between atoms and bosons. Thus, a character of order or 
disorder can be assigned to each system. 

In other words, a well specified PODI categorization emerges naturally. Specifically, the 
inequality a > (3 implies a character of disorder^ while ii a < (3 the system has the character 
of order and there is the possibility of gradual change of their character from order to 
disorder realized in the chain Atoms- Atomic Clusters- Nuclei-Bosons (from left to right), or 
equivalently from top to bottom in Table [H It is seen in Fig|6] that nuclei and bosons lie 
very close to the a-axis (disorder), atoms on the /3-axis (order), while for atomic clusters 
the disorder character is much stronger than that of order. 

The proper values of (a,/3) are obtained by an overall fitting of Ta,f}{N) with two free 
parameters to a single curve C{N) (no free parameters). Hence our statement that a quan- 
tum system is ordered or disordered comes from an overall comparison of the corresponding 
values of a and /3, i.e. for a range of values of A^. This enables us to find the optimal pair 
(a,/3). If one wishes to obtain values (a,/3) imposing the condition Vo,^f^{N) = C{N), for a 
fixed value of A^, then one has to fit a two parameter formula ra,i3{N) to just a single value 
C(A^), obtaining a dependence /5(a). We also observe that for all systems studied in this 
work the functional of disequilibrium D{N) is very similar to ra,j3{N) and C{N). The same 



observation was made in 



Iq . Thus a new measure of complexity emerges. 



The usefulness of our approach may be validated pragmatically, by future applications 
of PODI to several other quantum and classical systems. First, it should give reasonable 
results appealing to intuition. Second and more important, it is interesting to investigate 
cases where, instead of varying just the number of particles A^, one might study the effect 
on complexity of other relevant quantities of the systems as well. For example, the next step 
may be to study explicitly the effect of correlations between particles and/or the influence 
of various theoretical models of the systems. 

It is noted that so far there is no single or perfect quantitative deflnition of complexity. 
It is quite natural that any effort to provide a speciflc deflnition will be met with criticism 
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27, 



28 



29|. Comments on the validity of SDL and LMC measures can be found in Sec. 4 



of IJ], where welcome properties of any definition of complexity are described in detail. 

An answer to the question whether complexity shows an increasing or decreasing trend 
with A^, is dictated by the behavior of C{N). It is seen, that, at least for the models employed 
in the present paper, the trend of complexity is clearly decreasing for atomic clusters, nuclei 
and bosons, while it is increasing for atoms. An additional merit of the present work is 
the proper choice of the strength of disorder parameter a and the strength of order jS, 
by imposing the approximate condition T^^isiN) = C{N) for all values of A^. This is in 
accordance with a relativistic treatment of atoms, where C (N) shows an increasing behavior 



as well 



30 



311]. However, in our opinion, the latter result is not conclusive, because the 



corresponding studies are taking into account only the density distribution in position-space 
p(r). A proper, complete treatment of information entropy 5* and disequilibrium D should 
take into account both position- and momentum-spaces, indicated by relations S = Sr + Sk 
and D = Dr ■ Dk, employed in our present calculations. 

So far, {a,P) indices have been found for quantum systems. The same PODI procedure 
can be applied to quantify order versus disorder in classical systems, provided that they can 
be treated probabilistically. Such an application might be the use of probability densities in 
chaotic mappings, enabling researchers to quantify and observe the change of strengths of 
disorder a and order j3, while non-linear systems evolve through various phases of chaotic 
behavior and routes to chaos. Finally, one may, in principle associate to any probability 
distribution depending at least on one or more parameters describing a system, a PODI 
pair and a corresponding complexity character, under certain conditions. This case can be 
examined in a future work. 

The present study focuses on the dependence of complexity on just one parameter, i.e. 
the number of particles. The same method might be extended to assess the effect of other 
relevant quantities, which are of interest in the context of a specific level of description of 
a system. Another point, related to the latter one, is that in order to answer a specific 
question about organized complexity and the interplay of order versus disorder, one might 
need a sophisticated model of the dynamics of the system, or, in some cases, a simplified 
one, describing only salient but necessary features. This might be sufficient for a special aim. 
From this point of view, the merit of a model describing a system will be assessed, not only 
by its ability to reproduce exactly experimental results or simulate its behavior accurately. 
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but by being robust enough to uncertainties or errors in description. This is facihtated by 
the fact that an inequahty between a, (5 can be satisfied for a wider range of values, instead 
of an approximate equaUty a ~ /3. 

In brief, we propose a recipe to find a Pair of Order-Disorder Indices (PODI) in 

order to classify the order versus disorder contribution to complexity of any system (in our 
case quantum) described probabilistically. This enables us to assign to systems the category 
of order (a < /3) or disorder (a > /3) and classify them accordingly. 

We conclude that the atom is the only system out of four, studied here, which is ordered 
and growing in complexity with the number of particles (electrons). We conjecture that this 
is an information-theoretic explanation, supporting the fact that atoms arc used by nature, 
as suitable components to construct larger structures i.e. molecules and macromolecules. It 
is also a first step, with encouraging results, to fulfill quantitatively Wheeler's plan as the 
fundamental theory of the Cosmos. 
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FIG. 1: Atoms 
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FIG. 2: Atomic Clusters 




FIG. 3: Nuclei 




FIG. 4: Bosons 



17 



20 
15 
a 10 








20 
15 
alO 











5 10 

(a) Atoms 



10 



(c) Nuclei 



15 20 




15 20 



2 
15 
a 10' 












(b) Atomic Clusters 



20 
15 
a 10 







5 



10 



(d) Bosons 



10 15 20 
P 




15 20 



FIG. 5: Complexity regions in the a — (3 plane, showing a decreasing trend of rQ,,/3(iV) 
(white), a convex (grey) and an increasing one (black), for (a) atoms, (b) atomic clusters, 
(c) atomic nuclei, and (d) bosons. The approximate linear boundaries of the three regions 

are given in Table [B 
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FIG. 6: PODI values («,/?) for atoms, atomic clusters, nuclei and bosons. The a-axis 
corresponds to "disorder" while the /3-axis to "order". 
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